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We have grown single crystals of Na cc Ca ?/ Co02 and determined their superstructures as a function 
of composition using neutron and x-ray diffraction. Inclusion of Ca 2+ stabilises a single superstruc¬ 
ture across a wide range of temperatures and concentrations. The superstructure in the Na + layers is 
based on arrays of divacancy clusters with Ca 2+ ions occupying the central site, and it has an ideal 
concentration Na^rCai/rCoCb. Previous measurements of the thermoelectric properties on this 
system are discussed in light of this superstructure. Na^Cai/rCoCb corresponds to the maximum 
in thermoelectric performance of this system. 

I. INTRODUCTION 

Sodium cobaltate (Na cc Co02) has emerged as a ma¬ 
terial of exceptional scientific and technological interest 
due to its potential for thermoelectric applications that 
convert waste heat into electricity or allow solid state 
refrigeration 1 4 . In Na a ,Co02, the removal of sodium 
ions creates holes in the cobalt layers. The Seebeck co¬ 
efficient and power factor both increase with x, as the 
concentration of holes decreased 4 . Above a concentration 
of x ~ 0.85, however, phase coexistence occurs with the 
insulating x — l phase, causing a reduction in the power 
factor. The substitution of a divalent ion for sodium de¬ 
creases by one the number of holes in the cobalt layer. 

This offers the possibility to lower the hole concentra¬ 
tion and enhance further the thermoelectric performance. 

There are already promising results on polycrystalline 
samples doped with Ca and Sr, with an increase in the 
power factor over the comparable pure compouncP^. 

Neutron diffraction studies have demonstrated the cen¬ 
tral role played by superstructure on the thermoelectric 
properties of Na^ CoOj^. The structure of Na cc Co02 
comprises layers of conducting cobalt-oxygen octahedra 
separated by layers of sodium ions. The sodium ions 
can sit on one of two inequivalent sites; Na 2b sites di¬ 
rectly between successive cobalt ions, or displaced from 
these position on Na 2d sites. The 2d sites are pref¬ 
erentially occupied due to a small energy cost for sit¬ 
ting directly above cobalt arising from short-range repul¬ 
sion. Vacancies in the sodium layer have been found to 
form multi-vacancy clusters at certain fractional concen¬ 
trations. These clusters form by promoting sodium ions 
onto the 2b sites, creating rattling cages that can reduce 
thermal conduction through the material 10 . Coulomb 
interactions cause the vacancy clusters to arrange them¬ 
selves into long range patterns, causing rings of superlat¬ 
tice peaks in diffraction images. 

The electrostatic potential generated by these super¬ 


lattices templates the Coulomb landscape in the cobalt- 
oxygen layer, generating confined electronic pathways in 
the layers 8 . The combination of electron confinement 
and rattling cage-ions creates an ideal environment for 
the thermoelectric effect. 

Here we present the results of an x-ray and neu¬ 
tron diffraction study of superstructures in Ca-doped 
Na a ,Co02. We show that the Ca stabilises a single super¬ 
structure over a wide range of temperature and compo¬ 
sition, and we use the results to gain an understanding 
of how Ca-doping affects the thermoelectric properties. 

II. METHODS 

Single crystal boules of various compositions of 
Na ;E Ca ?/ Co 02 were grown using the floating zone 
technique 11 12 . Small, high quality crystallites were 
cleaved from the boules and studied using x-ray diffrac¬ 
tion with a molybdenum source, Xcalibur diffractometer 
(Agilent instruments) at Royal Holloway. A Cryojet5 
(Oxford Instruments) was used for temperature control 
and precise temperatures were calibrated by replacing the 
sample with a thermistor of similar size. A large, single 
grain of Nao. 7 Cao.iCo 02 was cleaved from its boule and 
studied using neutron diffraction on SXEP^at ISIS. 

The hexagonal reflections from the diffraction experi¬ 
ments were integrated using standard instrument meth¬ 
ods, and properties of the average structure were refined 
using the program Jana2006 14 . A reverse Monte Carlo 
(RMC) program was written to determine the vacancy 
ordering within the sodium layer and to establish the 
position of calcium ions. RMC uses simulated annealing 
to minimise the chi-squared difference between the cal¬ 
culated intensities arising from a given set of atomic po¬ 
sitions in the unit cell, and the experimental data. The 
positions of Ca and Na ions are changed with respect 
to the Co lattice, using the Metropolis algorithm in the 
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Canonical ensemble. Small displacements of all atoms, 
with respect to their original positions, are also allowed. 
When the fictitious temperature is slowly decreased, the 
configuration corresponding to the absolute minimum of 
this chi-square difference can be reached. Further details 
on this method can be found in Refs. [ 15 l 16 ] . 


III. DIFFRACTION EXPERIMENTS 

Three compositions of the calcium doped sys¬ 
tem were grown using the floating zone tech¬ 
nique: Nao.65Cao.iCo02, Nao.57Cao.i4Co02 and 

Nao.55Cao.2Co02- The samples were grown as large 
4-6 cm boules which were cleaved to find small, high 
quality crystallites. These were screened with the x-ray 
diffractometer and high quality samples were measured 
with longer exposures and greater coverage at several 
temperatures. Figure [l] shows reciprocal space cuts 
in the (ft, ft, 0) plane from a single-crystal sample of 
each composition. All three samples exhibit rings of 12 
satellite reflections around the principal hexagonal Bragg 
peaks. These additional peaks lay on the vertices of a 
hexagonal grid with spacing a*/7 and can be indexed 
using a supercell with lattice vectors: 

a' = 2a — b 
b' = a + 3b 


These supercell vectors index one of the two possible su¬ 
perlattice domains, where the second domain can be gen¬ 
erated by a reflection in either the {1,0,0} or {1,1,0}. 
The intensities of the superlattice reflections indicate 
equal populations of these domains. No broadening of the 
peaks was detected within instrumental resolution and, 
therefore, the superstructures were found to be ordered 
long range. The positions of these superlattice peaks are 
not consistent with any of the diffraction patterns ob¬ 
served previously in pure Na x Co 02 , signifying that a new 
superstructure has been found. All three samples have 
additional reflections halfway between the largest Bragg 
reflections due to contamination by A/2, in addition to 
this the Nao. 55 Cao. 2 Co 02 sample displays further reflec¬ 
tions between the primary satellites. These additional 
reflections lay at the vertices of an a*/4 hexagonal grid 
and can be indexed by a supercell with lattice vectors: 

a' = 2a 
b' = a + 2b 



FIG. 1. (Colour online) Reciprocal space cuts through the 
(ft, ft, 0) layer for different compositions of Na^Ca^CoCb, mea¬ 
sured by single-crystal x-ray diffraction, (a) Nao. 65 Cao.iCo 02 
shows a 12-fold ring of satellite peaks around the princi¬ 
pal Bragg reflections which can be indexed on a hexagonal 
grid with spacing a*/7. (b) Nao.57Cao.i4Co02 displays the 
same superlattice structure, but with stronger intensity, (c) 
Nao. 55 Cao. 2 Co 02 also displays the a */7 superlattice peaks, 
plus additional peaks that lay on a a*/4 lattice. 


These propagation vectors were first observed for 
Nao.sCoOp^, implying a coexistence between the 
Nao.5Co02 superstructure and the new a*/7 superstruc¬ 
ture. 

The intensities of the satellite peaks vary significantly 
for different / planes as shown in Fig. [2] Satellites occur 
primarily in 12-fold rings around the {1,0,/} reflections, 
with the modulation of peak intensity within the rings 
varying between even and odd planes. At high /, the peak 
intensity is diminished due to the x-ray form factor, how¬ 
ever peaks can still be observed and show that satellites 
around the {1,1,/} reflections gain intensity, especially 
in the (ft, ft, 11) plane. The increased satellite intensity 
and reduced intensity in the principal reflections in the 
(ft, ft, 11) plane is particularly striking and it is related to 
the fact that the c-lattice parameter is approximately 11 
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times the distance between the cobalt and oxygen planes. 

The diffraction patterns were measured as a function of 
temperature between 90 and 460 K. It was found that the 
satellite peaks do not change below room temperature, 
though they do eventually disappear at high tempera¬ 
ture, as illustrated in Fig. [3j The disappearance of the 
satellite peaks is attributed to the sodium layers becom¬ 
ing disordered. We did not observe any variation of the 
peak widths in this temperature range within our exper¬ 
imental resolution. The ordering temperature increases 
as the calcium concentration increases, and the order pa¬ 
rameter indicates a sharper transition for the y = 0.2 
sample. The a*/4 peaks did not vary at all within this 
temperature range. 

The peak intensities in Fig. [3] are normalised against 
the intensity of a nearby principal hexagonal reflec¬ 
tion in each case. The Nao.65Cao.iCo02 intensity is 
much weaker than the other compositions and it has 
a lower disorder temperature, with a fitted order pa¬ 
rameter giving Td = 379 =b 6K. There is little differ¬ 
ence between the initial intensity or the disorder temper¬ 
ature of Nao.57Cao.i4Co02 and Nao. 55 Cao. 2 Co 02 , with 
Td = 401 ± 2K and Td = 410 db 10 K respectively. This 
suggests that the a*/7 satellites in both samples origi¬ 
nate from a phase with a common calcium composition. 
The additional calcium in the Nao.55Cao.2Co02 sample 
is, therefore, likely to be part of the coexisting phase. 

A boule of Nao.7Cao.iCo02 was screened on SXD. A 
large single grain with intense, sharp principal Bragg 
peaks was cleaved from the boule and longer exposures 
of several hours were taken at 350K and 40K. It was not 
possible to grow large enough single grains for neutron 
diffraction studies of the superstructures with higher cal¬ 
cium concentrations. 

Superlattice peaks can be distinguished in the data by 
their intensities having hexagonal symmetry. Satellite 
peaks were observed with the same propagation vectors 
used to index the a*/7 superlattice observed with x-rays, 
indicating that this phase forms in the bulk material, 
as shown in Fig. [4] Several shorter exposures were per¬ 
formed at increasing temperatures and the disorder tran¬ 
sition observed was consistent with the Nao. 65 Cao.iCo 02 
measurement performed with x-rays, see Fig. [3j The 
superlattice pattern did not change as a function of tem¬ 
perature, indicating a single phase below the ordering 
transition temperature. 

In addition to the satellites, further peaks were ob¬ 
served close to the {110} principal hexagonal reflections. 
The d-spacing of these peaks was ~ 3.70 A, which 
matches the {220} reflection of calcium oxide (CaO). 
By comparing the integrated intensities for several re¬ 
flections with calculated structure factors, and averag¬ 
ing over domains, the relative proportion of CaO in 
Nao. 7 Cao.iCo 02 was determined to be 0.8% with a sta¬ 
tistical error of 0.2%. For this system to produce a hexag¬ 
onal pattern of points as found in these results, the cu¬ 
bic [111] direction of the CaO crystal would have to be 
parallel to the hexagonal [001] of Nao. 7 Cao.iCo 02 , and 


the cubic [110] direction of CaO would be parallel to the 
hexagonal [100] of Nao.7Cao.iCo02- In this direction the 
CaO structure can be described as successive hexagonal 
layers of calcium and oxygen ions with an ABC stacking 
sequence. The hexagonal environment of the main phase 
allows the CaO impurity to grow epitaxially within the 
Nao. 7 Cao.iCo 02 boule. 


IV. SINGLE CRYSTAL REFINEMENTS 

The principal hexagonal reflections in the multi¬ 
composition x-ray datasets were integrated using the 
three-dimensional (3D) profile method. The results 
of the refinements, as recorded in Table [T| show that 
the sodium occupancies are consistently lower than the 
nominal values. The calcium occupancies are much 
closer to the nominal values however, which is consis¬ 
tent with loss of the more volatile sodium ions during 
growth. For Nao.57Cao.i4Co02 the fitted Ca occupancy 
agrees with the nominal composition for the proposed 
Na 4 / 7 Ca 1 /7Co02 superstructure. This is consistent with 
the lack of any observed coexisting phases in this sample. 

The thermal parameters were refined as isotropic 
spheres. The data in Tableware consistent with increased 
static disorder away from the optimum composition. The 
ADPs for the sodium ions on the 2d sites are comparable 
to those obtained for the pure systeirPH, whereas those of 
the divalent calcium ions on the 2b sites are significantly 
lower. A full discussion of ADPs in terms of the lattice 
vibrations is beyond the scope of the present paper. 

The R-factors for these refinements show that there is 
excellent agreement between the refined models and the 
experimental data. Other refinements were attempted 
that changed the assumption that Ca ions occupy 2b sites 
and Na ions sit on 2d sites. These resulted in unrealistic 
concentrations of calcium or negative thermal parame¬ 
ters, and these models were, therefore, rejected. 


V. REVERSE MONTE CARLO SIMULATIONS 

Intensities from the superlattice reflections from 
Nao. 57 Cao.i 4 Co 02 were obtained from the room tempera¬ 
ture x-ray data by indexing a single superlattice domain 
and excluding the principal hexagonal reflections. The 
superlattice reflections arise purely from the superstruc¬ 
ture, whereas the principal reflections are only slightly 
modified, and they are also affected by other factors such 
as multiple scattering, or coexisting domains and phases. 
The propagation vectors used to index the rings of su¬ 
perlattice reflections converts to a real space supercell of 
seven hexagonal sodium cobaltate unit cells, and the con¬ 
centrations from the refinements meant that four sodium 
ions and one calcium were added to each layer. RMC cal¬ 
culations were performed using this set of intensities and 
supercell. Each iteration of the RMC program involves 
discrete hopping movements of the sodium and calcium 
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FIG. 2. (Colour online) (a) Cuts through the reciprocal space volume for Nao.57Cao.i4Co02 at integer l values, illustrating 
the variation of intensity in the superlattice peaks with l. Gaps in the data are a result of missing coverage, (b) Calculated 
intensities from the model generated by reverse Monte-Carlo, for comparison with (a). 


Occupancy Ui SO 

Target Composition a (A) c (A) Na Ca Na Ca Co O R w 

Nao.esCao.iCoOs 2.8398(6) 10.746(2) 0.56(1) 0.12(1) 0.019(2) 0.013(6) 0.0066(5) 0.006(1) 7.84% 

Nao.57Ca 0 .i4Co0 2 2.8403(4) 10.756(1) 0.529(5) 0.138(5) 0.010(1) 0.005(2) 0.0049(3) 0.0058(6) 3.96% 

Nao.55Cao. 2 Co02 2.8362(4) 10.761(2) 0.516(7) 0.172(9) 0.013(1) 0.014(3) 0.0062(4) 0.0050(7) 7.39% 


TABLE I. Refinement of Na ;r Ca 2/ Co0 2 principal Bragg reflections from x-ray diffraction data at 300K. The same model was 
used in each case: Space group: P^/mmc, Positions: Co (0,0,0), O (|, |, 0.0930(8)), Na (§, |,0.25), Ca (0,0,0.25). Co and 
O occupancies were fixed at full occupancy. 


ions to vacancy sites within the plane, as well as contin¬ 
uous movements of any ion by small increments in any 
direction. Simulated annealing of the system helped to 
avoid a false minimum solution, initially allowing ions to 
hop freely then slowly decreasing the fictitious tempera¬ 
ture and reducing the probability of allowing movements 
that would produce a higher % 2 . 

The RMC calculations led to the robust solution illus¬ 
trated in Fig. [5j with the calcium ion on a 2b position 
surrounded by sodium ions sitting on 2d sites. In the 
second sodium/calcium layer the same pattern emerges, 
however, the calcium ion is translated away from the loca¬ 
tion of calcium on the first layer. This translation relates 
to a significant reduction in % 2 . The fit for this model 


had a x 2 per degree of freedom of 5.82 and an R w value of 
12.9%. Figure [2|b) shows the calculated x-ray patterns 
for this solution indicating good agreement with the ex¬ 
perimental data in Fig. [2|a) across all l planes, including 
the variation of the satellite peak intensity between even 
and odd planes. The resulting structure has the space 
group P2i/m with unique axis c, atomic positions are 
given in Table \H\ 

The RMC solution also exhibits ordering in the distor¬ 
tions of the cobalt plane, with cobalt ions moving away 
from the closest calcium (which is sitting directly above 
the cobalt). This form of distortion was predicted using 
a Coulomb model for the pure systenP. The magnitude 
of the distortion is smaller for the system doped with 











































































5 



Temperature (K) 

FIG. 3. (Colour online) Temperature dependence of the su¬ 
perlattice peaks for three compositions of Na x Ca y Co02. Each 
peak intensity is normalised by a close principal hexagonal 
reflection, therefore indicating the relative intensity variation 
between compositions. Fitting an order parameter to each 
dataset gives disorder temperatures of Td = 379 zb 6K for 
Nao.65Cao.iCo02, Td = 401 d= 2 K for Nao.57Cao.i4Co02 and 
Td = 410 =b 10 K for Nao.55Cao.2Co02. The neutron data 
points indicate measurements during heating and cooling, 
showing no sign of hysteresis behaviour and are consistent 
with the x-ray measurement. 


calcium since the decreased hole concentration reduces 
the electrostatic potential for the multi-vacancy clusters. 
Figure [5ja) depicts the calcium doped divacancy struc¬ 
ture with cobalt distortions emphasized for clarity. It 
was not possible to resolve the variation in the Co-0 
bond lengths using our x-ray data and, therefore, we at¬ 
tempted to use the neutron data. 

RMC simulations using the neutron diffraction data 
were less sensitive to the superlattice, mainly due to a 
lower signal-to-background. The same divacancy super¬ 
structure was obtained, but the smaller scattering con¬ 
trast between the calcium and sodium ions made it diffi¬ 
cult to locate the positions of the calcium ions. 

For Nao.55Cao.2Co02 this divacancy phase coexists 
with the phase observed for Nao.5Co02 as shown in Fig. 
[6j In fact, the Nao. 5 Co 02 superstructure can also be 
viewed as stripes of divacancy cluster^. In this case 
it was not possible to perform RMC simulations due to 
the presence of contamination by higher order diffraction 
peaks. 


VI. DISCUSSION 

Analysis of the superlattice observed in all samples of 
Na^Ca^CoC^ indicates the presence of the superstruc¬ 
ture with an ideal concentration of Na^Ca^CoC^. 
For pure sodium cobaltate the corresponding superstruc¬ 



Wavevector Q _L b* (A 1 ) 


FIG. 4. (Colour online) A cut in the (h, k, 0) plane for 
Nao.7Cao.iCo02 SXD data at 40K. Satellite peaks appear in 
the same positions as found in the x-ray data, indicating that 
this is a bulk phase. Black circles and squares distinguish the 
two domains for the superlattice, and the red diamonds mark 
the {220} reflections of the epitaxial CaO impurity. Other 
non-indexed peaks are a result detector anomalies and noise 
resulting at low detector angles. 


Name 

X 

y 

z 

Site 

Col 

0 

0 

0 

2 a 

Co2 

0.1453(2) 

0.7156(2) 

0.00453(18) 

4f 

Co3 

0.2850(2) 

0.4295(2) • 

-0.00029(17) 

4f 

Co4 

0.4266(2) 

0.1421(2) 

0.00147(17) 

4f 

01 

0.0462(11) 

0.2425(12) 

0.0926(9) 

4f 

02 

0.2436(11) 

0.1951(11) 

0.5877(9) 

4f 

03 

0.1948(11) 

0.9557(12) 

0.0977(9) 

4f 

04 

0.3759(11) 

0.9031(11) 

0.5922(9) 

4f 

05 

0.3284(11) 

0.6601(12) 

0.0955(8) 

4f 

06 

0.5271(11) 

0.6205(11) 

0.6018(9) 

4f 

07 

0.8966(11) 

0.5176(11) 

0.0853(9) 

4f 

Nal 

0.2334(14) 

0.1921(15) 

0.25 

2 e 

Na2 

0.0592(14) 

0.2410(15) 

0.75 

2 e 

Na3 

0.3799(14) 

0.8966(15) 

0.25 

2 e 

Na4 

0.5226(14) 

0.6177(15) 

0.25 

2 e 

Cal 

0.1486(8) 

0.7151(8) 

0.75 

2 e 


TABLE II. Atomic coordinates of the Nao.57Cao.i4Co02 
structure, determined by RMC. The structure has the mon¬ 
oclinic space group P2i/m with unique axis c. The lat¬ 
tice parameters are a' = 7.5177(6)A, b' = 7.5201(6)A, c — 
10.7564(9)A, 7 = 119.986(8)°. 

ture has a sodium concentration x = 5/7. According 
to calculations of the ground state energies of multi¬ 
vacancy structures using a long-range Coulombic model 
with short-range repulsion, this is the only stable diva¬ 
cancy cluster in this composition rangd^. It is likely that 
the substitution of the divalent calcium at the centre of 
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FIG. 5. (Colour online) Results from the RMC calculation 
of Nao.5rCao.i4Co0 2 x-ray data, (a) The resultant structure, 
where blue balls are sodium ions, purple are calcium ions, the 
other balls are cobalt ions, where dark green is a large neg¬ 
ative distortion and light green is a large positive distortion. 
Oxygen ions have been omitted for clarity, (b) Annealing of 
X 2 during the calculation. Inset: Final comparison of calcu¬ 
lated (J c ) and experimental (J 0 ) intensities, where the black 
line displays I c — I 0 - 


the divacancy cluster has a stabilising effect since it low¬ 
ers the charge of the divacancy cluster. The fact that 
this superstructure does not change with temperature 
and forms over a range of sample compositions suggests 
that it is a particularly stable superstructure. 

Previous studies involving the doping of sodium cobal- 
tate with calcium have concentrated on the variation of 
transport properties relevant to the thermoelectric per¬ 
formance of the material and had not considered the role 
of the superstructure. This has led to a number of dif¬ 
ferent doping regimes being attempted though none have 
used exactly the concentration of the superstructure de¬ 
termined here. Kawata et al found that both the electri¬ 
cal resistivity, p, and the Seebeck coefficient, 5, increase 
with increasing calcium concentration 5 . Ono et al. found 
that the increase in both p and S was linear, so that over¬ 
all the power factor, P = S 2 /p, also increased with cal¬ 
cium concentration^. However, the power factor reaches 
a maximum value for Nao. 7 oCao.oi 75 Co 02 and above this 
calcium concentration, P decreases. In another com- 
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FIG. 6. (Colour online) Nao.55Cao.2Co02 can be described 
by a coexistence of the Nao. 57 Cao.i 4 Co 02 supercell in (b) 
and the x — 0.5 supercell in (c). The calculation of 
the (/i, k, 0 ) layer in (a) gives good agreement with the x- 
ray data from Fig. [ljc). Superlattice peaks generated by 
the Nao. 57 Cao.i 4 Co 02 supercell are indicated with circles, 
whereas the peaks produced by the x m 0.5 supercell are 
marked by triangles. 


position series, Li et al. found that the power factor 
first increased, and then decreased for Nao. 64 Cao.i 6 Co 02 
and higher calcium concentration^. In these reports the 
power factor is peaked very close to the optimum con¬ 
centration for the ideal Na^Cai^CoC^ superstructure. 
The Seebeck coefficient increases up to this composition. 
Above this composition, where we find coexistence with 
an additional insulating phase, the electrical conductivity 
decreases. 

For a full understanding of the thermoelectric perfor¬ 
mance of this material, the effect of the divacancy su¬ 
perstructure on the thermal transport must also be un¬ 
derstood, as this is a key component for calculating the 
thermoelectric figure-of-merit, ZT = TS 2 jpK , where k 
is the thermal conductivity. This contribution is beyond 
the scope of the current work but will be addressed in a 
future publication. 


VII. CONCLUSIONS 

In summary, both x-ray and neutron diffraction of 
calcium-doped sodium cobaltate have revealed a previ¬ 
ously unobserved superlattice that can be explained by 
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a superstructure with seven unit cells. An RMC pro¬ 
gram was used to find the ordering in both sodium 
layers as well as distortions in the cobalt plane in a 
model-independent manner, reliably generating a diva¬ 
cancy structure with calcium ions at the centre on a 2b 
site, with an ideal concentration Na 4 / 7 Ca 1 / 7 Co02. The 
same superstructure is observed across a broad tempera¬ 
ture range and in all compositions studied, though mov¬ 
ing away from the optimal concentration for this struc¬ 
ture results in weaker superlattice intensities and the 
emergence of additional phases. Comparison with ther¬ 
moelectric measurements in the literature indicates that 


this ideal concentration coincides with the maximum in 
the thermoelectric power factor for the system. 
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